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Abstract 

We study the branching ratio, CP-violating asymmetry, forward-backward asymmetry 
and the CP-violating asymmetry in the forward-backward asymmetry for the exclusive 
decay B Kt^t^ in the two Higgs doublet model with tree level flavor changing neutral 
currents (model III). We analyse the dependencies of these quantities on the neutral 
Higgs boson contributions and the CP parameter sin6 in the model HI. We observe that 
to determine the neutral Higgs boson effects, the measurements of the forward-backward 
asymmetry and the CP- violating asymmetry in the forward-backward asymmetry for the 
decay B Kt~^t~ are promising. 
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1 Introduction 



Rare B meson decays, induced by flavor changing neutral current (FCNC) b ^ s transition 
are the most promising research areas to test the Standard model (SM). Since these decays are 
induced at loop level in the SM, a comprehensive information can be obtained for the more 
precise determination of the fundamental parameters, such as Cabbibo-Kobayashi-Maskawa 
(CKM) matrix elements, leptonic decay constants, etc. Further, they shed light on the physics 
beyond the SM, such as two Higgs Doublet model (2HDM), Minimal Supersymmetric extension 
of the SM (MSSM) 0, etc. 

Among the rare B-decays, the ones which the SM predicts large branching ratio (Br) become 
attractive since they are measurable in the near future, in the existing and forthcoming B- 
factories. The B — > Kl'^l' , {I = e, /i, r) decay, induced by 6 ^ sl^l^ transition at the quark 
level, is one of the candidate. In the literature, there are various experimental studies on this 
decay 0-0. The 90% C.L. upper limits of Br{B — > Kl^l ) {I = e, fi) have been obtained at 
the order of the magnitude as 10~^, close the SM predictions. 



In g -pO[, this transition has been investigated extensively in the SM, 2HDM . In these 
studies, the neutral Higgs boson (NHB) exchange diagrams are not taken into account since 
the lepton-lepton-Higgs vertices are proportional to the mass of the lepton underconsideration. 
However, for I = r case, the mass can not be neglected since it is comparable with the 
6-quark mass and NHB exchange diagrams give considerable contributions to the physical 



quantities of such channels. In [^, B Kt^t^ process is studied in the model II 2HDM and 
the NHB effects are taken into account. It is shown that the Br ratio of the process is enhanced 
for large tanjS values and the NHB contributions become considerable. Recently, this decay 
has been analysed and the forward backward asymmetry has been studied in the constrained 
minimal supersymmetric SM [^. 

The forward-backward asymmetry ApB and the CP- violating asymmetry Acp are the phys- 
ical quantities which provide information on the short distance contributions. ApB does not 
exist in the SM and also in the 2HDM without NHB effects. However, with the addition of 
these effects, ApB created and in the model II 2HDM, its numerical value increases with the 
increasing value of the vertex factor tanjS (see for details). The sources of Acp are the 
complex CKM matrix elements or the Yukawa couplings appear beyond the SM. Since the CP 
violation in the SM is negligible and no complex coupling exists in the model II (I) version of 
the 2HDM, one can go further and choose model HI version of the 2HDM to get a measurable 
Ac p. 
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In our work, we study the exclusive B Kt^t~ decay in the model III. Since the r lepton 
mass is comparable with 6-quark mass and the new Yukawa coupling ^^t-,-, coming from the 
vertices tt Hq or tt Aq can be large, we include the NHB diagrams and test the amount of 
their contributions. We calculate the Br of the process and observe that it is sensitive to the 
NHB effects. Second, we get non-zero Acp, at the order of the magnitude 10~^, since the 
Yukawa couplings in the model III can be taken as complex. Aps appears in the case that the 
NHB effects are non-zero and, therefore, we study its sensitivity to the Yukawa coupling ^^^^ 
and the mass ratio, of the neutral Higgs bosons Hq, Aq. Finally we calculate the CP- 

asymmetry in ApB {Acp{Afb)) and see that it is at the order of the magnitude 10~^. Similar 
to ApB, Acp{Afb) can exist if the NHB effects are non-zero and can be used for testing the 
contributions beyond the SM. 

Note that the theoretical analysis of exclusive decays is more complicated due to the hadronic 
form factors, which brings an uncertainity in the calculations. The calculation of the physical 
observables in the hadronic level needs non-perturbative methods to determine the matrix 
elements of the quark level effective Hamiltonian between the hadronic states. This problem 
has been studied in the framework of different approaches such as relativistic quark model 



by the light- front formalism [jT2|, chiral theory p3[, three point QCD sum rules method [p^ 



effective heavy quark theory and light cone QCD sum rules 

The paper is organized as follows: In Section 2, we calculate the Br, Acp, ApB and 
Acp{Apb) of the exclusive B Kt^t~ decay. Section 3 is devoted to the analysis of the 
dependencies of the physical quantities given above on the Yukawa coupling the ratio 



- and the CP parameter sinO. In appendix, we give a summary for the model HI and the 
calculation of the matrix element for the inclusive h sl^l~ decay in this model. Furthermore, 
we give the explicit forms of the form factors and the parametrizations used in the text. 

2 The exclusive B Kl^l~ decay 

The exclusive B — > Kl^l^ decay is induced by the inclusive h — > sl^l^ process which has been 
studied in the literature extensively. Recently the h — > sl^l~ decay has been handled with the 



addition the NHB effects in the framework of the general two Higgs doublet model In 
the appendix we give a summary of the model underconsideration and present the effective 
Hamiltonian which is used for the hadronic matrix elements. 

The calculation of the physical quantities like Br, ApB-, Acp, etc., need the matrix elements 
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{K |s7^(l + 75)6| B), {K \s{l + 75)6! B) and {K {sia^^q-'il + 75)6! B) and they read as [|T 
{K\s^,il-^,)b\B) = ipB+PK),r{q') + q,f-iq') , 



rriB + rriK 

{K \sil + j,)b\ B) = -L((ml-ml)r{q')+q'f-iq')) . (1) 

nib 

Here Pb and px are four momentum vectors of B and K mesons, q = Pb~Pk is the momentum 
transfer. Using these form factors, the matrix element of the B — > Kl^l~ decay can be written 
as: 

M = -^^VtbYZ I [ApK, + Bq,] hn + [CpK, + Dq^] I^ld + F^U + F^hd^ (2) 

where the functions A, B, C, D, Fi and F2 are given in Appendix B. Using eq.(^ and making 
the summation over final lepton polarizations, the double differential decay rate is calculated 
as: 

^'"'H^rJ^'"'^'' !^^^^^^'^'' + l\C\'ml (Ami + 16m^) + 2\F,\'ml s 

m^niiil — r ~ s) 

+ 8Re{C*D) m|mf (1 - r - s) + 2\Fi\^ml sv^ + z(4i?e(A*Fi) m| mi v) 

- yAm|.^(|Ar + |Cn)| (3) 

where z = cosO , 6 is the angle between the momentum of £ lepton and that of B meson in the 
center of mass frame of the lepton pair, A = 1 + — 2r — 2s — 2rs, i; = w 1 — ^ ~ 

V B B 

and s = In the light lepton case, namely I = e,u, the NHB effects are negligible and new 

TUb 

Wilson coefficients Cg^ and Cq^, appearing in the form factors Fi and F2, almost vanish (see 
27| for details). However, for r lepton, the NHB effects can give considerable contributions to 



dT 

dsdz 



the physical quantities -Br, ApB, Aqp, etc. Even if Aqp is possible without these effects, they 
play the main role in the existence of ApB- Further, Acp{Apb) exists when NHB effects are 
non-zero. Therefore, we concentrate on Acp{Apb) in addition to the quantities App and Aqp- 
Using the definitions 



we get 



and 



A„ = 5 , (6) 



Idsl 2A2A3\i v ml (1 - i)Im(XSM) 
^- = ^ Jl^ -' « 



where 



+ iA,FP+A,Fi'^)Re{^§^J + AM'^)Y (8) 

and Q is obtained from eq. @ by an integration over z and s of the terms in the curly bracket. 
A in eq. reads as 

A = vy/xi^mlX{Al + Al + C^){l-j) + 16ml{CD\^+{C\ + D^s)) 
+ 8m, F^^\VXC + 2sD) + As (jFap + v'^\Fi\'^) 

+ \^nm\"« XAl{l-j)+A (Ff )2 sv') + 2 Rei^^^J (A m| A, A,{1 AfI"^ ^ . 

+ Am^y/XC fP +8mrsDF^'^^)'^ (9) 

The exphcit forms of the functions Ai, C, D, f'-^\ i = 1,2, 3; j, /c = 1, 2 appearing in eqs. (H) 
and (|^) are given in Appendix B. 

Finally Acp{Afb) can be defined as 

AcAA.b) = 7 4" (10) 
where A^^^^ is given in eq. (j^) and ApB can be calculated by making the replacement — *• 

Notice that, during the calculations, we take into account only the second resonance for 
the LD effects coming from the reaction h — > sipi st^t~ , where i = 1, ..,6 and divide the 
integration region for s into two parts : < s < g^^^ {:m^^+Qm) ^ ^ q _ 7^)2 

where = 3.686 GeV is the mass of the second resonance 
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3 Discussion 



The model III induces many free parameters, such as where i,j are flavor indices and 
they should be restricted using the experimental results. Now, we would like to present the 
restrictions we use for these parameters, in our numerical calculations. Since the neutral Higgs 
bosons, ho and Aq, can give a large contribution to the coefficient Cj which is in contradiction 



with the CLEO data \28 



Br{B X,7) = (3.15 ± 0.35 ± 0.32) 10"^ , (11) 



the couplings C.^isi'^ = d,s,b) and db can be assumed as negligible to be able to reach the 



conditions (.^bb^Nis << 1 and (.^db^Nds << 1- (see the appendix of ||29| for details). Using 



^N,bb^N,is ^ '^'■'-'^ SAT.db SAT.ds 

also the constraints coming from AF = 2 mixing, the p parameter and the CLEO data we 
have: 

^N,tb ~ , ^^^ij ~ 0, ij = d,s quarks . 

This assumption permits us to neglect the contributions coming from primed Wilson coefficients 
which are related with the chirality flipped partners of the operator set (see for details) 
since the Yukawa vertices are combinations of ^§ sb ss- Finally, we only take into account 

the Yukawa couplings ^^^j^, ^N,bb ^^.^r- Notice that, for the coupling at first, we do 

not introduce any constraint and we will try to predict an upper limit by usin the present 
experimental measurements. At this stage we introduce a new parameter 6 with the expression 

FU t*D _ \FU F*D \ -id (^cy\ 
'iN,tt^N,bb — \^N,tt^N,bb\^ ■ y^"^) 

Here, it is possible to take both i^ii and or any one of them complex. In our work, we 
choose i^ii as real and ^^j,^ as complex, namely ^^^^ = iC^bfel The phase angle 6 leads to a 
substantial enhancement in neutron electric dipole moment and the experimental upper limit 
on neutron electric dipole moment dn < 10~^^e-cm thus places a upper bound on the couplings: 
-Im{^^ < 1.0 for mH± ^ 200 GeV ^|. 



mtnif, 



In this section, we study the dependencies of the Br, CP asymmetry Acp, forward-backward 
asymmetry ApB and CP asymmetry in forward-backward asymmetry Acp{Afb) of the decay 



In our analysis we restrict IC^'^'^^I in the region 0.257 < \Cj^^ < 0.439, coming from CLEO 



B Kt^t on the selected parameters of the model III {$,ntt^ phase angle 



measurement (see ||32| for details). With this restriction, an allowed region for the parameters 
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^Ntti ^Nbb ^1 is found. Throughout the numerical calculations, we respect this restriction, 
the constraint for the angle 9 due to the experimental upper limit of neutron electric dipole 
moment and take \rtb\ = \ jn^\ < 1, the neutral Higgs mass rriHo = 100 GeV^, charged Higgs 
mass mH± = 400 GeV, the scale ^ = nib- Here, we also give the input values used in the 
calculations, in Table (|^). 



Parameter 


Value 




1.40 (GeV) 


rrib 


4.80 (GeV) 


rrir 


1.78 (GeV) 


em 


129 


At 


0.04 




3.96 ■ 10-13 (GeV) 




5.28 (GeV) 


mt 


175 (GeV) 


mw 


80.26 (GeV) 


mz 


91.19 (GeV) 


J^QCD 


0.214 (GeV) 


as{mz) 


0.117 


sinOy/ 


0.2325 



Table 1: The values of the input parameters used in the numerical calculations. 

In figs. ID (H) we plot the Br of the decay B Kt^t^ with respect to the Yukawa 
coupling ^^^^^ (the ratio ^) for ^^ ^.^ = 40 mj,, rriho = 70 GeV, rriAo = 80 GeV {^n,tt = 10 m^, 
m^Ao = 80 GeV). The Br is restricted in the region between solid lines (dashed lines) for 
Cj^^ > {Cj^^ < 0). Fig. |1] shows that B r is sensitive to the NHB effects especially for 
case. For increasing values of ^^t-^, Br can take even two orders of magnitude larger 
values compared to the case where no NHB is taken into account. From this figure, it is possible 
to predict the upper limit of the coupling C,n,tt^ ^n,tt < 20 GeV, respecting the experimental 
upper limit, Br{B Kl^l^) < 0.5 x 10"^ , (/ = e,/i) with the assumption that the 
Br{B — s> Kt^t^) is not so much different than the previous one. For Cj^^ < the possible 
numerical values lie near 10"'' and the Br is not sensitive to the NHB effects. Note that, the 
Br in the SM is 1.0610"^ and in the model HI without NHB effects are between upper and 
lower limits (1.05-1.08) 10-^ ((0.95-0.97) 10-^) for C^^^ > {Cj^^ < 0). As shown in Fig. g 
the Br is also sensitive to the ratio for C^-^-^ > case and it increases for the larger values 
of the ratio. Furtermore, the experimental value of the Br{B Kl^l^) shows that the mass 
values of the neutral Higgs bosons and Aq are not far. 



6 



Fig. |] is devoted to the sin9 dependence of Acp including NHB effects. Here, Acp is 
restricted in the region between sohd hnes (dashed hnes) for Cj^^ > {Cj-^-^ < 0). It is at the 
order of the magnitude 10~^ and increases with the increasing values of the parameter sin6 as it 
should be. For Cj^-^ > 0, the possible values of Acp have the same sign (here minus), however 
for C^^^ < it can vanish or it can have both signs. We also present the sinO dependence 
of Acp without NHB effects, in Fig. ^ For this case Aqp is greater as a magnitude and the 
restriction region is larger compared to the previous one, expecially for Cj^^ > 0. As it can be 
seen from eq. (^, the addition of NHB effects reduces the magnitude of the Acp since NHB 
contributions enter into expression in the denominator but not in the numerator. 

In Figs. 1^ (|[), we plot the App of the decay B Kt^t~ with respect to the Yukawa 
coupling e^^,, (the ratio ^) for = 40 m^, nih, = 70 GeV,mA, = SOGeV = 10 m,, 

m^o = 80 GeV ). Ap p is restricted in the region between solid lines (dashed lines) for Cj^'^ > 
{C'/f < 0). Since App appears only with the NHB effects, it is a good candidate for testing the 
existence of them. App is at the order of the magnitude 10~^ and increases with the increasing 
values of f^,, for C^-^-^ < as shown in Fig. |[ The behavior of App is different for Cj^^ > 
since it slightly decreases when ^^^.^ increases. In addition to this, App is sensitive the ratio 
^ for Cf^ > and increases as a magnitude with the decreasing ratio (see Fig. However 
for Cj^^ < App is not sensitive to the ratio Further, it has negative sign for both 

Cj" > and C^" < 0. 

Finally, we present the CP violating asymmetry in App {Acp{App)) in Figs. |^ and |] since 
this parameter exists when NHB effects are non-zero and can play an important role in the 
determination of those effects. Fig. |^ shows dependence of Acp{App) for Cj-^-^ > (lies 
between solid lines) and C^-^-^ < (lies between dashed lines). This quantity is not so much 
sensitive to and can be at the order of the mag nitude 10-^ for C^^^ < 0. It can have 

both signs or can vanish for this case. For C^^^ > 0, the numerical value of Acp{Apb) can 
increase up to 10~^. Here, the SM Higgs Hq part of the NHB effects causes to have large values 
for Acp{App) and the part which contains neutral Higgs bosons beyond enters into expression 
destructively. Fig. || represents sin9 dependence of Acp{App). As shown in this figure, the 
possible values of Acp{App) for Cy-^-^ > have the same sign and they are non-zero for nonzero 
sinO, however for C^^^ < Acp{App) can vanish or exist with both signs. 

Now we would like to present our conclusions. 

• The Br of the exclusive process B Kt^t^ is at the order of the magnitude 10^'' for 
the SM and model HI withouth the NHB effects. However, including the NHB effects and 
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taking large values of the neutral coupling it is possible to enhance the Br more 

than one orders of magnitude compared to the one calculated in the SM. 

• It would be possible to predict the upper limit of the coupling ^^t-^, ^^tt < 20GeV, 
respecting the experimental upper limit, Br{B —>■ Kl^l") < 0.5 x 10~^ , {I = e, /i) ||^. 

• Using the complex Yukawa coupling (^^^^ it is possible to get a CP violating asymmetry 
Aqp at the order of the magnitude 10"'^, which is a measurable quantity. With the 
addition of NHB effects the magnitude of Acp decreases. 

• ApB is another physical quantity which exists when the NHB effects are non-zero. The 
calculations show that ApB is at the order of the magnitude 10~^ and the experimental 
measurement of this quantity can give strong evidence about the existence of NHB effects 
and the physics beyond the SM. 

• Finally, the CP asymmetry in App can appear also with the NHB effects and it is another 
physical quantity which can be used for testing the existence of the NHB effects. We cal- 
culate this quantity at the order of the magnitude 10~^ and its experimental measurement 
can give important clues about physics beyond the SM 

Therefore, experimental investigations of these physical quantities ensure a crucial test for 
the new physics beyond the SM. 
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Appendix 



A The model III and the inclusive b sr^r decay 

In the SM and model I and II 2IIDM, the flavour changing neutral current at tree level is 
forbidden. However, such currents are permitted in the general 2HDM, so called model III 
and it brings new parameters, i.e. Yukawa couplings, into the theory. These couplings are 
responsible for the interaction of quarks and leptons with gauge bosons, namely, the Yukawa 
interaction and in this general case it reads as 



Cy = vfiQiL<PiUjR + v^,QiL<PiDjR + &ij2UjR + CSQi^^^ 



+ r]j^JiL<j)iEjR + ^^JiL(j)2EjR + h.c. 



(13) 



where L and R denote chiral projections L{R) = 1/2(1=1=75), 0^, for k = 1, 2, are the two scalar 
doublets, QiL {Ul) are quark (lepton) doublets, UjR and Djr are quark singlets, Ejr are lepton 



singlets, rjfj^, and are the matrices of the Yukawa couplings. The Flavor changing (FC) 
part of the interaction is given by 



J^Y,FC = i^jQiLhUjR + &iL(t>2DjR + ^f^liLhEjR + h.c. 



With the choice of 0i and 02 



1 

71 







and the vacuum expectation values. 



< 01 > = 



, 02 



J_( V2H+ 



; < 02 >= 



(14) 



(15) 



(16) 



the SM and beyond can be decoupled. In eq. QI^) the couplings ^^'^ for the FC charged 
interactions are 



^ch ~ ^neutral VcKM , 
^di — ^CKM ineutral , 

where ^neutral Q is defined by the expression 

cU,D _ (-irU,D\-lcU,D-,rU,D 



(17) 



^In all next discussion we denote ^^^trai 
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Here the charged couphngs appear as hnear combinations of neutral couplings multiplied by 
VcKM matrix elements (see ^0|] for details). 

Now, we would like to present the procedure to calculate the matrix element for the inclusive 

b —>■ ST^T^ decay briefly: 

• Integrating out the heavy degrees of freedom, namely t quark, W^, charged Higgs boson 
if^, and neutral Higgs bosons Hq, Hi, H2 bosons in the present case and obtaining the 
effective theory. Note that Hi and H2 are the same as the mass eigenstates ho and Aq in 
the model 111 respectively, due to the choice given by eq. (plSj) . 



• Taking into account the QCD corrections through matching the full theory with the 
effective low energy one at the high scale fi = my/ and evaluating the Wilson coefficients 
from my/ down to the lower scale /i ~ 0{mf,). 

In the 2HDM, neutral Higgs particles bring new contributions to the matrix element of the 
process b —>■ st^t^ (see ||2^) since they enter in the expressions with the mass of r lepton or 
related Yukawa coupling C,ntt- Besides, there exist additional operators which are the flipped 
chirality partners of the former ones in the model 111. However, the effects of the latter are 
negligible since the corresponding Wilson coefficients are small (see Discussion part). Therefore, 
the effective Hamiltonian relevant for the process b st^t~ is 

neff = -A^VttV*^^Y.c^{^^)o^{^^) + T.CQXf^)Q^if^)] > (i9) 

where Oi are current-current {i = 1,2), penguin {i = 3, ...,6), magnetic penguin {i = 7,8) and 
semileptonic {i = 9, 10) operators. Here, Ci{fi) are Wilson coefficients normalized at the scale 
/i. The additional operators Qi{i = 1, .., 10) are due to the NHB exchange diagrams and Cq^{^) 



are their Wilson coefficients (see |^ for the existing operators and the corresponding Wilson 
coefficients). 

B The form factors and the functions appearing in the 
expressions 

We parametrize the fuctions A, Fi and F2 as 

A = Ai + iA2 + ^§^bbA3 
p _ I ED p(2) 

F2 = + (20) 
10 



with 



A, 



2Re{Clf^)t- ^^^^^ C^.ff 



f° ^0 



(1) 



mB + mKmbTrit 3 
-12/23 (m| - m^)/+ + m%sf- 



rrib 



^ 1 

NM) 



(2) 



-12/23 



{m% - m\)f + + m%sf~ 



N,bb 



dx 



l-x 



^(2) 



-12/23 



[m%-m]^)f + + m%sf 



rrib 



-12/23 



dx 



dx 



l-x 



dyC^'M^bb)) , 
dy (C^mhtf) + c^'M'^, 



l-x 



dyC^l^iU)) , 



(21) 



^b^Nfib 

where the formfactors /+, and /t are calculated in the framework of the light cone QCD 



sum rules and represented in the pole forms as ||26 

0.29 



r 
f 



1 - 



23.7 

0.21 



24.3 
0.31 



:22) 



23 



The other functions B, C and D appearing in eq. (Q) are 



B 
C 
D 



niB + rriK 



2Cio/^ , 



(23) 
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Figure 1: Br as a function of for fixed ^^^^ = 40 mf,, ruhg = 70GeV, 112^0 = SOGeV and 

sin9 — 0. Here Br lies in the region bounded by solid (dashed) lines for C^^^ > {C^^^ < 0). 
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Figure 2: Br as a function of for fixed thaq — SOGeV, ^ 

sinO = 0. Here Br lies in the region bounded by solid (dashed) lines for C^^^ > {C^-^^ < 0). 
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Figure 3: Acp as a function of sin9 , for fixed ^^t-^ = 10 m,-, ^^bi, = 40 m^, m/j,, = 70GeV, 
and thao — ^^GeV. Here Acp lies in the region bounded by solid (dashed) lines for C^^-^ > 
{C/^ < 0). 
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